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Oscillating a 2 -dynamos and the reversal phenomenon of the global 
geodynamo 
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Abstract. A geodynamo-model based on an a-effect which has been computed under conditions suitable for the geody- 
namo is constructed. For a highly restricted class of radial a-profiles the linear a 2 -model exhibits oscillating solutions on a 
timescale given by the turbulent diffusion time. The basic properties of the periodic solutions are presented and the influence 
of the inner core size on the characteristics of the critical range that allows for oscillating solutions is shown. Reversals are 
interpreted as half of such an oscillation. They are rather seldom events because they can only occur if the a-profile exists 
long enough within the small critical range that allows for periodic solutions. Due to strong fluctuations on the convective 
timescale the probability of such a reversal is very small. Finally, a simple non-linear mean-field model with reasonable 
input parameters based on simulations of Giesecke et al. (2005) demonstrates the plausibility of the presented theory with a 
long-time series of a (geo-)dynamo reversal sequence. 
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1. Introduction 

Paleomagnetic measurements show that the Earth's magnetic 
field exists for more than 10 9 years with nearly the same mag- 
nitude (Kono & Tanaka 1995). The process that is responsible 
for the production of this field is called the geodynamo and 
essentially takes place in the fluid outer part of the Earth's 
core. The magnetic field is dominated by a dipole which - 
as the most characteristic feature - from time to time "starts 
to oscillate" and changes its polarity from one sign to the 
other. This phenomenon has been called reversals and typi- 
cally lasts 10 3 - 10 4 years (Bogue & Merill 1992). The dura- 
tion of such "oscillating" phase is extremely short compared 
with the time between consecutive reversals. In fact the aver- 
age time between two reversals is about 50 times longer than 
the duration of the reversal itself. The appearance of the re- 
versals seems to be chaotic rather than periodic (Krause & 
Schmidt 1988), and the probability of a reversal during a cer- 
tain time span can be described by a Poisson distribution. On 
geological time scales (~ 10 7 — 10 8 years) the average rate 
of reversals changes (Merill et al. 1996) and there exist even 
very long periods were no reversal occurs - so called super- 
chrons. 
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Closely related to the reversal phenomenon are so called 
excursions a kind of aborted reversals, where the polarity be- 
gins to change but, instead of executing a full transition, the 
dipole returns to the original polarity. Excursions occur about 
ten times more often than reversals. 

Deviations from a perfectly axisymmetric field become 
manifest in the tilt of the dipole axis with respect to the ro- 
tation axis (currently 11°) and in the non-axisymmetric field 
components in terms of localized flux patches. Such field pat- 
terns in average exhibit a common directed drift motion, the 
so called westward drift (see e.g. Bloxham & Jackson 1989; 
Bloxham & Jackson 1992). Roughly simplifying the value of 
this westward drift amounts approximately 0.3°/year. 

Simulations of the 3D MHD-equations (that describe a 
thermal/chemical driven turbulent flow of a conducting fluid 
in the Earth's outer core and the magnetic field that is in- 
duced by this flow) have been able to reproduce some of the 
observed features of the Earth's magnetic field (Glatzmaier 
& Roberts 1996; Kageyama & Sato 1997; Christensen et al. 
1998; Kuang & Bloxham 1999). Unfortunately such calcula- 
tions are computational very expensive. The time periods that 
can be covered are rather short compared to the time scales 
on which for example changes of the mean reversal rate oc- 
cur. In order to examine geodynamo-models in matters of the 



2 THE EQUATIONS 



statistics of the reversal phenomenon mean-field models re- 
main indispensable. 

A further unsolved issue is the influence of the small 
scales. Global MHD simulations are restricted in the achiev- 
able parameter regime and in the affordable spatial resolu- 
tion. These limits prevent from resolving the actual scales 
of the turbulence in the fluid outer core, and for reasons of 
numerical stability unphysical large values for the viscous 
losses have to be adopted. Therefore the smaller scales are 
artificially damped, and properties and influence of the small 
scale turbulence remains unsure. 

Sarson & Jones (1999) developed a 2.5D model to exam- 
ine the reversal mechanism in detail. Their general picture 
is a "large scale" ail-dynamo mechanism, where a strong 
zonal flow and meridional circulations are responsible for the 
dynamo action. Reversals are induced by fluctuations of the 
meridional flow. The authors reclaim that a opposite dipole 
polarity can evolve if the velocity remains long enough in a 
regime that allows for oscillating magnetic fields. 

Statistical properties in a one-dimensional afi-dynamo 
model where reversals are triggered by a fluctuating a-effect 
have been examined by Hoyng et al. (2002). Their model was 
able to predict some basic features of the Sint-800 data and 
although there have been some contradictions in general their 
mean-field approach seems to be selfconsistent. 

Other models, based on the absence of differential rota- 
tion (or any shear flows) in the Earth's core, interpreted the 
geodynamo as an a 2 -dynamo (Steenbeck & Krause 1966). It 
is known that the spherical a 2 -dynamo "almost always" pos- 
sesses stationary axisymmetric magnetic field solutions (for 
scalar a-effect) or longitudinally drifting non-axisymmetric 
solutions (Rudiger et al. 2003) and therefore such models 
have difficulties to explain the reversal phenomenon. There 
are few exceptions of the rule that a 2 -dynamos with scalar 
a-effect are stationary. 

Fearn & Rahman (2004) solved the Navier-Stokes equa- 
tion and a mean-field induction equation for an a 2 -dynamo 
model with a radial dependence of the a-effect given by 
a oc sin7r(r — i?i n ). This radial profile leads to a vanish- 
ing a-effect at the boundaries of the fluid outer core but the 
a-effect does not show any zero within the interior. They ob- 
tained non-linear periodic solutions if the a was larger than 
a certain critical value. In contrast to the solutions of Sar- 
son & Jones (1999) their results were strongly influenced 
by the non-linear back-reaction of the Lorentz force on the 
flow which serves as a saturation mechanism for the mag- 
netic field. 

Without considering any mean flow, Stefani & Gerbeth 
(2003) found oscillating a 2 -dynamos in case that the a-effect 
(uniform in 8) changes its sign in radius. Already Soward 
(1974) with his quasi-linear approximation for rotating con- 
vection in layers with uniform density found that the a-effect 
strongly varies with depth: it is negative (positive) in the 
lower (upper) part of the convection layer - well described 
by a radial sinus-function. Giesecke et al. (2005) have shown 
that such profiles indeed result from numerical simulations 
of the convection in the outer fluid core where the density 
stratification is very small (see Sect. [3}. 



In the present paper we shall show that radial sin-profiles 
of the a-effect lead to oscillating a 2 -dynamos but already a 
slight deviation from this profile provides stationary modes. 
Combining the principle properties of the a-effect from the 
calculations of Giesecke et al. (2005) with a simple a 2 - 
dynamo, a mean-field model is constructed that exhibits ir- 
regular reversals induced by a fluctuating a-effect. 



2. The equations 



Taking the induction equation 
dB ( 

— - = Vx\uxB-ri\7xB 
dt V 



(1) 



and split magnetic field B and velocity u in a mean part, (B), 
(u) and a fluctuating component B', u' the mean magnetic 
field (B) = B - B' is determined by 



9(B) 
dt 



V x («) x (B) +£ -77V x (B) 



(2) 



with £ = (u' x Bj as the mean electromotive force (EMF) 
and rj the (molecular) magnetic diffusivity. The components 
of the EMF are usually given by 



Ei = aij(Bj) + p i]k dk{Bj). 



(3) 



The tensor a^ correlates the EMF due to turbulent motions 
with the large-scale magnetic field, including the effects of 
anisotropy. In the simplest case the tensor /3jjfc is related to 
the turbulent diffusivity by ftijk = r\T^ijk which is the case 
that we shall discuss here. In the following all mean flows 
(u) are neglected and we end up with a mean-field induction 
equation that describes an a 2 -dynamo: 



d(B) 
dt 



V x a(B) - ?7tV x (B) 



(4) 



where ijt fj is implied. Equation l|4} together with a pre- 
scribed a-effect that depends on the radius r and the latitude 
angle 8 is solved numerically using an explicit finite differ- 
ence scheme in two dimensions in spherical coordinates. The 
standard resolution is 64 x 64 grid points in both radial and 
latitudinal directions. A perfect conductor is assumed to exist 
at the inner core boundary which is justified because of the 
absence of turbulent motions in the solid inner core so that 
the diffusivity is significant smaller inside the solid core than 
inside the fluid outer core. However, a finitely conducting in- 
ner core affects the behavior of the magnetic field as it has 
been shown e.g. by Hollerbach & Jones (1993). At the outer 
boundary a vacuum is simulated by increasing the magnetic 
dissipation by a factor of 10. The details of the numerical 
realizations of these boundary conditions are described by 
Rudiger et al. (2003). To estimate general properties of the 
linear mean-field model at first a quenching mechanism is 
abandoned to avoid the complicated questions that are asso- 
ciated with the non-linearities. Only for the long time simula- 
tions in Sect.|4]an equilibration mechanism is used to prevent 
the field from growing to infinity. 



3.2 Oscillating a 2 -dynamos 



3 GEODYNAMO a -EFFECT 



3. Geodynamo a-effect 
3.1. General properties 

Fig-fflshows a typical radial dependence of the a-effect com- 
puted from local simulations of weak stratified and fast ro- 
tating magnetoconvection by Giesecke et al. (2005). The a- 
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Fig. 1. Radial dependence of the a-effect for A = - — — — = 

2/ipilrj 

1 . The dashed thick line denotes the time average and the grey 
lines show fluctuations in time. 

effect on the northern hemisphere shows a maximum (mini- 
mum) close to the upper (lower) boundary and the cross-over 
takes place exactly in the middle of the convective unsta- 
ble layer. The radial profile is almost perfectly antisymmetric 
with respect to the middle of the layer so that the net a-effect 
(integrated over r) approximately vanishes. 

The dotted lines divide the domain in three parts. In the 
outer zones the exact determination of the a-effect is difficult 
because close to the boundaries strong gradients of the mag- 
netic field require the exact knowledge of the turbulent dif- 
fusivity ?7t for an calculation of the a-effect from Eq. (0. In 
the central part - between r — 0.48 and r = 0.87 - the field 
gradients are negligible and thus the result should be more 
reliable. But even if the presented profile comes with some 
uncertainties, the qualitative overall behavior can quite well 
be described by a ~ sin r, where the argument of the sin 
must be chosen in a way that the a-effect disappears at the 
inner and the outer core boundary. If a density stratification 
is included in the simulations, the a-effect cross-over moves 
more and more towards the bottom of the box (see Riidiger 
& Hollerbach 2004 their Fig. 4.23). 

3.2. Oscillating a 2 -dynamos 

The general properties of the above presented a-effect are 
used as an input for a global axisymmetric a 2 -dynamo. The 
model is a spherical shell with the inner radius Ri n and the 
outer radius i? ou t. For the present date Earth the radius of the 
solid inner core is given by i?; n = 1222 km and the radius 
of the fluid outer core is given by i? ou t = 3480 km. In all 
simulations the radius of the outer core is scaled to 1 and 
to maintain the correct ratio, the radius of the inner core is 



scaled to 0.35. In the following the a-tensor is antisymmetric 
with respect to the equator (~ cos 9). The "standard profile" 
of the a-effect is given by 

a(r, 9) = ao cos 9 sin 27r— ^— I . (5) 

V -Rout — -Rin / 

This equation is slightly modified to vary the amplitude in the 

upper (lower) half of the sphere and the zero-crossing of the 

a-effect. We start with a strict sin-profile in radius of the a- 

effect, i.e. with a cross-over i?o in the middle between i?j n = 

0.35 and i? ou t = 1 and equal amplitudes. If supercritical the 

field grows exponentially and the resulting dynamo oscillates. 

The condition to Rq for the excitation of the oscillation 

of a 2 -dynamos is very strict. The oscillations only exist for 

Ro « 0.67. ..0.70. If the radial a-profile does not lie in the 

narrow area indicated in Fig. |2] (top) then the dynamo does 

not oscillate. The periodic solution due to the strict radial sin- 
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Fig. 2. Critical domains for the radial profile of the a-effect 
that lead to oscillating solutions. Top: critical interval for the 
location of the zero. Bottom: critical interval for the magni- 
tude of the lower amplitude 

profile of the a-effect can also be suppressed by a variation of 
one of the amplitudes. Numerical experiments were made by 
multiplication of the lower part of the sin-profile with a factor 
A. Oscillating solutions are found with 8 A ~ 0.1 (see bottom 
of Fig. |2ji . Of course, a similar result would hold for the upper 
part of the radial sin-profile but the general property of the 
model is now known: The oscillating solution only exists if 
the deviations of the actual a-effect profile from the profile 
given by Eq. l|5} are very small. 
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3.4 Oscillation period and critical dynamo number 



Due to fluctuations a simple radial sin-profile of the ot- 
effect is rather seldom. A possible oscillation of the dynamo 
only happens if the profile of the a-effect lasts long enough 
within the critical range that allows for periodical solutions. 
The minimum time which must be covered by the (critical) 
radial a-profile in order to excite (half of) an oscillation has 
been estimated from the simulations and is given approxi- 
mately by 



imin ~ 0.3 • T di ff 

with the diffusion time Tdiff defined by 



Tdiff 



v 



(6) 



(7) 



If this minimum time is not reached by the radial profile of 
the a-effect the stationary solution would not start to become 
oscillating. 

In all simulations the time the dipole needs for the transi- 
tion from one polarity to the other is of the order of Tdiff ■ If we 
assume -q ~ 2 • 10 4 cm 2 /s (for molten iron under conditions 
in the fluid outer core) we retrieve Tdiff ~ 10 5 years. To keep 
the model consistent with the observed duration of a reversal 
(10 4 years) we have to assume a turbulence-induced enhance- 
ment of the magnetic diffusivity by one order of magnitude: 

VT « 10 • rj. (8) 

With i]t ~ 20m 2 /s the diffusion time reduces to 10 4 yrs, the 
typical duration of a reversal. The ratio r/T/r) ~ 10 seems 
not to be totally devious and was also the result of a rough 
estimation of Giesecke et al. (2005). A comparable value has 
been presented by Hoyng et al. (2002) who independently 
determined a turbulent diffusion time of 10000 — 15000 years 
from the analyzation of autocorrelation functions from the 
Sint-800 observations. However, it should be kept in mind, 
that this are very rough estimations. 

3.3. Field pattern of a reversal 

Figure [5] shows the temporal behavior of the magnetic field 
projected on a meridional plane during one oscillation. The 
left-hand side of each panel shows the isolines of the toroidal 
field component where solid (dashed) lines denote field direc- 
tions clockwise (counterclockwise). The arrows on the right 
hand side represent the poloidal field component. The length 
of the arrows is scaled with the field length. The magnetic 
field is concentrated near the rotation axis which is a con- 
sequence of the latitudinal cos 6>-dependence of the a-effect. 
The solution is of dipolar parity. Regarding the toroidal com- 
ponent in the northern hemisphere the considered reversal 
cycle starts close to the rotation axis with a clockwise ori- 
ented toroidal field in the upper half of the outer core and 
a field of opposite sign starts in the lower part of the outer 
core. The outer part weakens (1,2) and is replaced by a grow- 
ing toroidal field of opposite sign (2,3,4). Between the two 
belts of equally counterclockwise (clockwise) oriented com- 
ponents that determine the appearance of the magnetic field 
in the northern (southern) hemisphere in the middle of the re- 
versal sequence (5,6) a new opposite directed field appears 
(7) and pushes away the lower toroidal component (8). The 



reversal is completed in panel (9) when the lower part of the 
sphere is completely filled out by this new reversed oriented 
field. Note that the poloidal field component already shows 
the reversed polarity in part (6) of Fig.[3]and only undergoes 
minor changes in the remainder of the sequence. 



3.4. Oscillation period and critical dynamo number 

Figure|4]shows the oscillation period (dashed line, in units of 
Tdiff) and the critical dynamo number (solid line) 



£fcrit 



-tin 



(9) 



in dependence of the location of the zero of the a-effect {Rq). 
£fcnt determines the minimum amplitude of the a-effect at 
which dynamo-action occurs. The vertical dotted lines indi- 
cate the transition between oscillating and stationary solu- 
tions. The critical dynamo number is always larger for the 
oscillating solutions. This is not a surprising result. The os- 
cillating solutions have smaller scales than the stationary so- 
lutions so that the Ohmic losses are larger for the oscillat- 
ing solutions. The immediate consequence is that for given 
amplitude of a the oscillating solutions are less nonlinear 
than the non-oscillating solutions. In the center of the crit- 



o 





Q- 



no oscillations 




no oscillations . 




1 
1 

1 1 


- \ 




1 i 
\ 1 

! t 











o ES 

o 



0.68 

Rn 



Fig. 4. The characteristics of the linear a 2 -dynamo model in 
dependence of the location of the cross-over point Rq in the 
radial a-profile. 



ical interval the oscillation-period is nearly constant for Rq 
and increases strongly if the zero of the a-effect is close to 
the upper or lower boundary that separates the oscillating so- 
lutions from the stationary states. The values of the critical 
dynamo number for the symmetric and antisymmetric (with 
respect to the equator) axisymmetric modes (SO, AO) and for 
the first non-axisymmetric modes (Al, SI) are specified in 
Tabled Since the dominant part of the Earth's magnetic field 
is a dipole the AO mode is of profound interest. Note that the 
AO-mode (leads to a dipole solution) and the SO-mode (leads 
to a quadrupole solution) coincide within our numerical ac- 
curacy. 



3.5 Influence of the inner core size 
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Fig. 3. The magnetic field configuration of the cyclic solution with a strict radial sin-profile of the a-effect. The time in units 
of the diffusion time is printed in the lower right corner of each snapshot. 



Table 1. Critical dynamo number for different basic dynamo 
modes for a taken from Eq. (Rq = 0.675) 



AO 


Al 


SO 


SI 


C*= rit 17.04 


17.67 


17.04 


17.67 



3.5. Influence of the inner core size 

The solid inner core of the Earth is growing on geological 
time scales (hundreds of millions of years). Due to the spe- 
cific thermodynamic conditions in the Earth's interior the liq- 
uid iron first freezes out at the center of the Earth because the 
melting point of iron decreases faster with increasing pres- 



sure than the temperature increases towards the center. The 
appearance of a solid inner core few billion years ago resulted 
in important changes of the physical conditions and processes 
that dominate the flow in the fluid core. It is obvious that the 
size of the solid inner core also should have a strong influence 
on the oscillations of the a 2 -model that has been presented 
above. We restrict our examination to the geometric effects 
that arise from different sizes of the inner core and do not 
consider the changes in the turbulence that are associated e.g. 
with the emerging compositional convection. Figure[3]shows 
the critical profiles of the a-effect that lead to oscillating so- 
lutions for different ratios of R ou t/Rin- For increasing size 
of the inner core - indicated by the solid vertical line on the 
left side - the critical interval becomes smaller (indicated by 
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the dashed line). At the same time the center of this interval 
moves closer to the center of the fluid outer core (denoted 
by the dotted line). For Ri n = 0.1 (bottom curve) the criti- 
cal cross-over-points of the a-profiles are clearly located in 
the upper half of the fluid outer core, whereas for i?j n = 0.7 
(top curve) the critical profile is nearly a perfect sin-profile as 
given by Eq. l|5} where the zero is located exactly in the mid- 
dle of the fluid outer core and the width of the critical interval 
has become very small. 
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Fig. 5. Critical a-profiles for oscillating solutions for dif- 
ferent sizes of the inner core. From top to bottom: i?j n = 
0.7,0.5,0.35,0.2,0.1. The dashed lines confine the critical 
interval and the dotted line denotes the center of the fluid 
outer core. 

Although the critical interval becomes smaller with in- 
creasing size of the inner core - making it more difficult to 
excite an oscillation - one can speculate that the overall prob- 
ability for a reversal in case of a fluctuating a-profile might 
be approximately constant (or at least > 0) for all sizes since 
the center of the interval moves towards the middle of the 
fluid outer core which is the preferred location for the zero 
cross-over of the a-effect (see Fig.^. 

3.6. Non-axisymmetric modes 

The a-effect is a nontrivial tensor if the rotation is fast: 
ct-zz —* for f2 — > oo (in cylindrical ccordinates, see Moffatt 
1970; Riidiger 1978; Busse & Miin 1979) where a zz refers to 
the a-effect in cylindrical coordinates and z denotes the di- 
rection parallel to the rotation axis. The tensorial structure 
of the a-effect is now taken into account, i.e. relation the 
a zz — is used in the dynamo equation. 

Figure [6] shows C™* (solid line) and the drift period 
(dashed line) for the lowest non-axisymmetric modes. Again 
the results for the symmetric mode (SI) and the antisymmet- 
ric mode (Al) coincide in both quantities. The dotted vertical 
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Fig. 6. C° rlt of the non-axisymmetric modes for the case 
a zz = 0. The corresponding drift periods are given by the 
dashed line. 

lines denote the critical interval for oscillating solutions for 
the AO- respective SO-mode for the (scalar) isotropic a-tensor 
as described in the previous section. 

For the zero of the a-effect below Rq 0.705 we obtain 
a westward drift. Above Ro 0.71 both modes show an east- 
ward directed drift motion. The characteristic drift time scale 
is approximately (0.20.. .0.30) -T^ig for the westward drifting 
modes and (0.40...0.50)-Tdiff for the eastward drifting modes. 
Note that the ratio of the drift timescale (0.2rdiff) to reversal 
timescale (raiff) obtained from the simulations coincides with 
the same ratio observed for the Earth where the timescale of 
the westward drift is ~ 2000 years and the duration of a rever- 
sal ~ 10000 years. However, these timescales are estimations 
with a large uncertainty (e.g. data for the reversal time reach 
from 100 to several 10000 years). 

The critical dynamo numbers for the axisymmetric modes 
(SO, AO) are larger than C"'* of the Al/Sl-modes so that 
the solution would be dominated by these non-axisymmetric 
modes. The axisymmetric modes with the higher eigenval- 
ues are oscillating. We know that this constellation is in con- 
tradiction to the observations of the Earth's magnetic field 
which is dominated by a stationary dipole-part. The non- 
axisymmetric modes only lead to smaller contributions that 
are manifested in the dipole tilt and the drifting field pat- 
terns. The interaction of the nonaxisymmetric modes (for 
anisotropic a) and the oscillating modes (for isotropic a with 
cross-overs) is still an open question. 

4. Irregular reversals induced by a fluctuating 

a-effect 

In the following the complications that appear from the non- 
axisymmetric modes are ignored. The effects of the full non- 
trivial a-tensor will be treated in a subsequent paper. 

It is known from numerous calculations that the a- 
coefficients are rather noisy quantities (e.g. Ossendrijver et 
al. 2001; Giesecke et al. 2005). Figure0shows the maximum 
(minimum) of the a-effect in dependence of the time in units 
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of the turnover time or advective timescale taken from the 
simulations of Giesecke et al. (2005) 

Rout ^Mn 



7"adv — 



(10) 



(where v! is the turbulent rms-velocity). 




Fig. 7. Fluctuations of the maximum and the minimum of the 
a-effect in dependence of time (in units of the turnover-time 
r a dv) taken from the results shown in Fig.[0 



The amplitude of the a-effect in the outer part of the shell 
is slightly larger than in the inner part and the strength of the 
fluctuations amounts approximately 10% of the average. It is 
obvious that the timescale of the fluctuations is given by r a dv 
If we assume a typical value for the Earth's fluid outer core: 
u' k, 5 • 10 _4 m/s the timescale of the fluctuations r ac i v is 
related to the diffusive timescale Tdis by 

Tadv « (0.01. ..0.02)- T difI (11) 

where Tdis is estimated by Eq 0. Thus the duration which 
is necessary for the a-effect to stay in a critical state - given 
by Eq. - is about 15... 30 times longer than the time scale 
on which the a-effect fluctuates. This indicates that a reversal 
clearly must be a very seldom event. It is also typical for the 
presented theory of the reversal phenomenon of the geody- 
namo that practically never a realization of the a-profile may 
exist so long that a complete oscillation can happen. 

For the long time calculations we adopt an isotropic a- 
effect given by Eq. (|5J and add fluctuations for the magnitude 
and for the location of the zero. For simplicity we assume 
equal averages for the upper and lower amplitude which both 
vary independently. The fluctuations are described by a Gaus- 
sian distribution with a standard deviation a of 10% of the 
average. The fluctuations of the zero cross-over are slightly 
larger. We adopt an average value of Rq = 0.675 (the middle 
of the fluid layer) with a standard deviation of a = ±0.100. 
According to Eq. il lb the actual values for the fluctuating 
quantities are updated each O.OlTdig . 

The dynamo number is given by C a = 20. This is clearly 
overcritical and therefore the non-linearities introduced by a 
local quenching function for the a-effect given by 

a ( B ) = a o— ^2 (12) 
1 + r> 



might also have some influence on the solution. Indeed test 
calculations show that increasing C a - corresponding to a 
stronger driven and thus more non-linear dynamo - reduces 
the probability of a reversal. 

Figure [8] shows the time dependence of the radial mag- 
netic field at some point in the spherical shell for a long time 
calculations that spans 2000 diffusion times (with the actual 
time scaling this corresponds to 20 million years). The mag- 
netic field reverses irregularly. Both polarity states occur with 
nearly the same probability and field strength. The distribu- 
tion of time-periods between each reversal can be described 
by an exponential function e _At /*, thus the reversals are in- 
dependent randomly occurring events (see Fig.|9jl. The num- 
ber of short life-times between consecutive reversals is over- 
estimated because it is difficult to distinguish between excur- 
sions and reversals. Within the presented 2D-simulations the 
only possibility to characterize a reversal is a sign change of 
a field component. Other properties, e.g. an slight increase of 
the tilt of the dipole axis followed by an immediate decrease 
back to the original state as an indication for an excursion, 
are intrinsically not available. Here we filtered out all events 
where the field changes its sign only "slightly" and for a very 
short time (< 0.5 Tdig). This means that an original dipole- 
state recovers very fast after the magnetic field just touches 
the zero in Fig. [8] 

For the chosen set of parameters the mean reversal rate 
is about a factor of 5 higher than the rate observed for 
the present date Earth. Short-time test-calculations also have 
been performed with slightly different values for t/t /r/. This 
is the crucial relation that specifies the a-effect fluctuation 
timescale with respect to the turbulent diffusion time. The 
choice of this timescale fixes the time span for which the a- 
effect remains in a certain state and therefore determines the 
probability of a reversal. Decreasing ?7t/?7 increases the (tur- 
bulent) diffusion time which results in a decrease of the re- 
lation T a dv/idiff (assuming that r a d v is fixed for the Earth's 
case). Therefore the timescale of the fluctuations of the a- 
effect is shorter with respect to the diffusive timescale and 
the probability of a reversal decreases. The results show that 
already a slightly reduced value for T a d v /Tdiff significantly 
increases the mean time between consecutive reversals. 
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Fig. 9. Number of realizations of life-times of a single polar- 
ity state. 
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Fig. 8. Reversals of the magnetic field. 173 events have been identified as reversals which leads to a mean polarity life-time 
of approximately 1 1 diffusion times. 



5. Discussion 

We have shown that an a 2 -dynamo is able to exhibit periodic 
solutions where the dipole polarity changes on a timescale of 
the order of the diffusion time Tdiff. A crucial contribution for 
the existence of an oscillating a 2 -dynamo is an a-effect that 
strongly depends on the radial coordinate r. The most impor- 
tant property of such a-effect is a sign change within the field 
producing domain and a certain symmetry between positive 
and negative magnitude. Only a highly restricted class of (ra- 
dial) a-profiles leads to oscillating solutions. It turned out that 
the radial profiles of the a-effect that have been calculated by 
Giesecke et al. (2005) under conditions that are suitable for 
the Earth's fluid core provide the essential characteristics that 
are necessary to construct such oscillating dynamos. 

Long time simulations over 20 million years retrieve solu- 
tions which show irregular reversals. The model is very sen- 
sitive to the choice of parameters that describe the statistical 
properties. In principle it should easily be possible to adopt 
values for the fluctuating quantities (C a , a and T ac i v /Tdiff) in 
a way that would result in a time series which exhibit the ba- 
sic properties of the geodynamo. The most important charac- 
teristics are the mean reversal rate, the distribution of polarity 
states, the distribution of the mean time between reversals and 
the ratio between number of excursions to number of rever- 
sals. An exact adjustment of the model at this stage is beyond 
the purpose of this work because further effects that affect the 
behavior of the field as well have not been considered. We ne- 
glected the influence of a finitely conducting inner core. Test 
calculations have shown that this does not prevent the solu- 
tions from oscillating, only some properties like the positions 
of the critical intervals might be slightly changed. Possible 
changes in the behavior of the convection (changes in ampli- 
tude and strength of the fluctuations) have been ignored just 
as the (possible) existence of large-scale flows. 

An unsolved issue is the inclusion of the non-trivial com- 
ponents of the a-effect. Currently it is not possible to re- 
trieve the characteristics of the non-axisymmetric parts of the 
Earth's magnetic field. The Ansatz a zz = - based on the- 
oretical reasons - resulted in a dynamo that is dominated by 
the non-axisymmetric modes (and the axisymmetric modes, 



dipole and quadrupole are strongly suppressed). A resulting 
realistic solution should be described by a combination of dif- 
ferent modes as it seems to be the case for the geodynamo. 
Indeed, observations of the magnetic fields of other planets 
or moons in the Solar system show, that various manifes- 
tations of field configurations are possible. The nearly per- 
fect dipole field of Saturn or the non-axisymmetric dominated 
fields from Uranus and Neptune are some extraordinary ex- 
amples. 
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